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ABSTRACT 

Context. The Seyfert 1.8 galaxy ESO 113-GOlO had been observed for the first time above 2keV by XMM-Newton during a short 
exposure (~4ks) in May 2001. In addition to a significant soft X-ray excess, it showed one of the strongest (in EW) redshifted FeKa 
lines, at 5.4 keV. 

Aims. We present here a long (100 ks) XMM-Newton follow-up of this source performed in November 2005, in order to study over a 
longer time-scale its main X-ray properties. 

Methods. We use both timing analysis (Power Spectra Density analysis, rms spectra, flux-flux analysis) and spectral analysis which 
mainly focuses on the Fe Ko- line complex. 

Results. The source was found in a higher/softer time-averaged flux state, and timing analysis of this source reveals strong, rapid 
variability. The Power Spectral Density (PSD) analysis indicates (at 95% confidence level) a break at 3.7^] !J x 10""* Hz. This cut-off 
frequency is comparable to those measured in some other rapidly- variable Seyferts, such as MCG-6-30- 15 and NGC 405 1 . From the 
mass-luminosity-time-scale, we infer that Mbh ranges from 4 x lO^-lO' Mq and the source is accreting at or close to the Eddington 
rate (or even higher). The existing data cannot distinguish between spectral pivoting of the continuum and a two-component origin for 
the spectral softening, primarily because the data do not span a broad enough flux range. In the case of the two-component model, the 
fractional offsets measured in the flux-flux plots increase significantly toward higher energies (similar to what is observed in MCG- 
6-30-15) as expected if there exists a constant reflection component. Contrary to May 2001, no significant highly redshifted emission 
line is observed (which might be related to the source flux level), while two narrow emission lines at about 6.5 keV and 7keV are 
observed. The S/N is not high enough to establish if the lines are variable or constant. As already suggested by the 2001 observation, 
no significant constant narrow 6.4 keV Fe Ko- line (EW<32 eV) is observed, hence excluding any dominant emission from distant cold 
matter such as a torus in this Seyfert type 1.8 galaxy. 

Key words, galaxies: Seyfert - galaxies: active - X-rays: galaxies - accretion discs - quasars: individual: ESO 1 13-GOlO 



1. Introduction 

The Seyfert ESO113-G010 (z=0.0257) was observed for the 
first time in X-rays, in December 1995, by ROSAT, revealing 
an X-ray luminosity of about 2.5 x lO'^-'ergs"' (assuming 
Ho-75 km s "' Mpc~' and (7o=0 .5) and some indication of X-ray 
variability ("Pietsc h et alj Il998h . An optical follow-up (2.2 m 
ESO/MPG telescope at la Silla) performed in November 1996 
done by the same authors to characterize this bright galaxy 
led to the conclusion that this Seyfert is of type 1.8, and 
measured logL;r/R=44.34erg s"', log Lb =43 .97 erg s"\ and 
logLx,ROSAT=43.40erg s"' . This object was observed for the 
first time above 2keV with XMM-Newtorj on May 3, 2001 
during a short exposure of about 4 ks. That observation revealed 
a soft excess below 0.7 keV and, more interestingly, a narrow 
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emission line at 5.4 keV (rest-frame) , most probably origi nating 
from a redshifted iron FeKa line ( Porquet et al.l l2004b). The 
line was detected at 99% confidence, calculated via Monte 
Carlo simulations which fully account for the range of energies 
where a narrow iron line is likely to occur. Such naiTow spectral 
features in the 5-6 keV energy range were discovered with 
XM M-Newton and Cha ndra in sever al other AGN: e.g. NGC 
35 16 (iTumer et al.ll2002) NGC 73 14 dYaaoob etalllool. Mkn 
766 (lTurneretal.1 l2004h . AX J0447-0627 (iDella Ceca et al.1 
120051) . The redshifted FeKa Une discovered in ESO 113-GOlO 
had an EW of ~ 265 eV, making it one of the strongest (in 
EW) redshifted iron lines observed to date in an AGN. The 
energy of the redshifted line could indicate emission from 
relativistic (0. 17-0.23 c) ejected matter m oving away f r om th e 
observer, as proposed for Mrk 766 by iTurner et al.l (l2004l) . 
Alternatively, if the classification as type 1.8 object is correct, 
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then emission from a narrow annulus at the surface of the 
accretion disc is unlikely due to the very small inclination 
angle (i.e. less than 10°) required to explain the narrow, 
redshifted line in this intermediate Seyfert galaxy. However, 
emission from a small, localized hot-spot on the disc, occurring 
within a fraction of a complete disc orb i t, could also explai n 
the redshifted Hne (e.g.. lTurner et al.ll2002l:lDovciak et al.ll2004l) . 

Here, we present a lOOks XMM-Newton follow-up observa- 
tion of ESO 113-GOlO. In sectionlH we present the long XMM- 
Newton analysis performed in November 2005, and the corre- 
sponding data reduction used for this work. Section [3] presents 
the timing analysis: light curves, PSD analysis, rms-spectrum, 
and flux-flux analysis. Section |4] presents the spectral analysis, 
focusing on the study of the FeKa line complex between 6.4- 
7 keV. Finally, we summarize and discuss the main results of the 
present work (section|5]l. 

2. XMM-Newton observation and data reduction 

XMM-Newton observed ESO 113-GOlO starting on November 
10, 2005 (orbit 1085) with exposure times of ~103ks and 
~ 102ks for flie MOS , and PN respectively. The EPIC-MOS 
cameras (iTurner et al .! 2001) operat ed in the L arge Window 
mode, while the EPIC-PN camera ( Striider et al.l 1200 ID oper- 
ated in the standard Full Frame Window mode. Background flar- 
ing was negligible during most of the observation and there- 
fore the whole observation was used, except when mentioned 
in Section [3] The EPIC data were re-processed and cleaned 
using the XMM-Newton SAS version 6.5.0 (Science Analysis 
Software) package. The netresulting livetimes were 102. 1 ks and 
90.6 ks respectively for the MOS and PN cameras. Since pile-up 
was negligible. X-ray events corresponding to patterns 0-12 and 
0-4 events (single and double pixels) were selected for the MOS 
and PN, respectively. Only good X-ray events (with FLAG=0) 
were included. The low-energy cutoff was set to 300 and 200 eV 
for MOS and PN, respectively. The source data were extracted 
using a circular region centered on the source position of diame- 
ter of 50"Q, 25", and 35" for MOS 1, MOS2 and PN respectively. 
The diameter was smaller for MOS2 compared to MOS 1 due to 
a dead pixel column about 26" from the source in the MOS2 
data. ESO 113-GOlO was by far the brightest X-ray source in 
this 30' EPIC field-of-view. Background data were taken from 
box regions on the same CCD as the source (excluding X-ray 
point sources). The xspec vll.3 software package was used for 
spectral analysis of the background-subtracted spectrum using 
the response matrices and ancillary files derived from the SAS 
tasks RMFGEN and arfgen. 
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Fig. 1. Soft (S: 0.2-2 keV) and hai'd (H: 2-10 keV) 500-s resolu- 
tion (background-subtracted) PN light curves of ESO 1 13-GOlO 
(upper and middle panels, respectively), together with the hard- 
ness ratio between the two bands (lower panel). The grey stripes 
correspond to the Icr error bars. For comparison the count rate 
and hardness ratio level found during the 4ks May 2001 obser- 
vation are reported (dashed lines). 



energy range are 2.42+0.02cts s ' and 3.71+0.01, for the 
2001 and 2005 observations, respectively. The corresponding 
Hardness Ratio, HR (=(H-S)/(H-i-S)) is also displayed (Fig. □ 
lower panel). The error bars for the hardness ratio are obtained 
using Gaussian error propagation (see e.g., Pietsch et al. 2005|). 
The 2005 HR values show that the source was on average in a 
softer state compared to May 2001. There was no very strong 
variation amplitude of the HR, except within the first lOks (with 
<AHR>~0.07). 



3. Timing Analysis 

PN light curves (time bin of 500 s, and error bars of Icr) clearly 
show fairly rapid variability with relatively large amplitude, 
both in the soft (S: 0.2-2keV) and hard (H: 2-lOkeV) energy 
ranges (Fig. [T] upper and middle panels, respectively). The 
corresponding PN time-averaged count rate observed during 
May 2001 (duration of ~4ks) is reported (Fig. [1] dashed 
lines) and shows that the source was in a higher flux state for 
virtually all of the 2005 observation when compared to the 2001 
observation. The PN average count rates in the 0.2-12 keV 

' For the spectral analysis (§4) the extraction region radius was taken 
to be the same for MOSl and MOS2 (i.e. 25 "), in order to co-add both 
spectra. 



To quantify this variability better, we carried out timing anal- 
ysis using the power spectral density function (PSD), the rms- 
spectrum, and flux-flux analysis. To maximize the signal-to- 
noise we combined events from all three instruments, re-scaling 
the background regions used for the MOS detectors to match the 
source/background area ratio of the PN. We used only data from 
times between 8 ks and 95 ks from the start of the observation, 
to remove the times of non-negligible background flaring for se- 
cure timing analysis, and merged the good time intervals of the 
different instruments in order to use only data received simulta- 
neously from all three instruments. The combined events file was 
used to extract light curves in different bands to make the PSDs 
and rms-spectra. Although the instrument responses differ, in all 
cases (including rms-spectra) we normalize by the mean flux in 
each band, so the effect of the spectral response is taken out. 
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Fig. 2. 0.3-lOkeV PSD of ESO 113-GOlO, including the ratio of 
the data to the best-fitting power-law plus constant noise compo- 
nent model (model shown as solid line). 
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Fig. 3. 0.3-lOkeV PSD of ESO 113-GOlO, including the ratio 
of the data to an exponentially cut-off power-law (with low- 
frequency slope a - -1), plus constant Poisson noise compo- 
nent (model shown as solid line). 



3.1. PSD Analysis 

We first consider the PSD for the light curve in the 0.3-lOkeV 
energy range, which is shown in Fig. ID The PSD represents the 
average of the (squared) amplitude of variations as a function of 
temporal frequency. The data are then fitted with a simple power 
law plus the flat level of variability due to Poisson noise expected 
given the error bars. A simple power-law plus noise is not a very 
good fit (x^ - 24.1 for 1 1 d.o.f.), and is rejected at > 98% confi- 
dence (for power-law slope a = -2.1). Allowing the noise level 
to be free improves the fit to - 18.9 for 10 d.o.f., which 
is still rejected at > 95% confidence. Finally, leaving the noise 
level free while allowing a break in the PSD improves the fit: fix- 
ing the low-frequency slope to -1 (simila r to the low-frequenc y 
PSD shapes observed in other AGN, e.g. McHardv et al.ll2005h . 
we find a break frequency of (5.9 + 1.7) x 10"^ Hz (errors 
correspond to 90% confidence for a single interesting parame- 
ter), with an upper limit on high-frequency slope of a < -2.9. 
The corresponding = 11.9 for 9 d.o.f., which, according 
to the F-test, is a significant improvement over the unbroken 
power-law model at the 95% level. We also find a good fit 
(X^ = 13.9 for 10 d.o.f.) for an exponentially cut-off power- 
law, with low-frequency slope fixed to - 1 and cut-off frequency 
(3.7+} °) X 10-4 Hz, which we show in Fig.[3l 

We caution that the possible break or cut-off in the PSD is 
only marginally significant, but note that, in support of these 
models, the product of frequency and power, which is con- 
stant on the slope=-l part of the PSD is ~0.01, which is 
of the same order as that observed in other AGN. In com- 
parison with other AGN, the measured cut-off frequency is 
comparable to those measured in some other rapidly- variable 
Seyfe r ts, such as MCG-6-3 0-15 (8 x lO^^ Hz. iMcHardv et al.1 
200l IVaug han et ail l2003b ). and NGC 4051 (8 x lO"'* Hz, 
McHardy et al.l2004 . The correlation between AGN X-ray vari- 
ability PSD break time-scale and black hole mass is now well- 
established. It appears that AGN PSD break time-scales can be 
linearly scaled down by the black hole mass to match the breaks 
observed in the P SDs of stellar mass black holes in X-ray binar y 
systems (e.g. Ma rkowitz et al.1 120031: lUttlev & McHardvll2005l) . 
Recently, McHardy et al. (2006) showed that the residual scat- 
ter in the mass-time-scale relation could be explained if time- 
scales also scale with luminosity or equivalently, accretion rate. 
Combining AGN and stellar mass black hole data for Cyg X-1, 



McHardy et al. framed the best-fitting relation for break time- 
scale logTe - 2.11ogMBH-0.981ogLboi-2.33, with Tb inunits of 
days and Mbh and bolometric luminosity Lboi in units of 10^ M© 
and 10"^^ erg s"' respectively. The best-fitting PSD break time- 
scale we measure, 3.7 x 10 ^ Hz, corresponds to = 0.03 days. 
For a range in Lboi from lO^'*'^ - 10'*^ erg s"' (i.e. assuming a 
minimum equal to the sum of FIR, optical and X-ray luminosi- 
ties reported bv Pietsc h et al.|[r99 8', and in this work), w e infer 
from t he mass-luminosity-time-scale relation of McHar dv et al.l 
(l2006h that Mbh ranges from 4 x 10''-10^ Mq and the source is 
accreting at or close to the Eddington rate (assuming even higher 
Lboi implies that the AGN is a super-Eddington accretor). 

We also measured PSDs using smaller subsets of the ob- 
served energy range, but due to the relatively weaker S /N of 
these PSDs, the energy dependence of the PSDs could not be 
well constrained. 

3.2. The rms-spectrum 

The rms-spectrum quanti fies the spectral shape of the variable 
part of the emission (see IVauphan et alj l2003al for a review). 
Specifically, one measures the root-mean-square (rms) variabil- 
ity in light curves made from a number of narrow energy bands 
or 'channels', and then normalizes each rms value by the mean 
in that channel to obtain the fractional rms (also known as Fvar)- 
A plot of Fvar versus energy yields the shape of the variable part 
of the spectrum relative to the mean spectrum (i.e. independent 
of instrumental spect ral response). The error bars are determined 
using the formula of 'Vaugh an et al.l (12003 al) . We show the mea- 
sured rms-spectra in Fig.lH obtained using light curves binned to 
5 ks (top panel) and 500 s (bottom panel). The rms-spectra ob- 
tained with 500 s binning naturally includes faster fluctuations 
than for the 5 ks binned data, and the Fvar measurements are 
accordingly larger Fvar appears to decrease above about 2 keV 
for both bin timescales. Therefore the soft X-ray excess down 
to 0.3 keV appears more variable compared to the hard energy 
range. This is different from the typical rms-spectra consistent 
with a broad peak in the variability amplitude between 0.8 and 
2ke V, and a m uch lower variability amplitude below .8 keV 
(e.g., iFabian et al. 2002; Inoue & Matsumoto 2003; Pont i et al.1 
120061) . In addition, there is a hint of a dip in the 500 s bin 
rms-spectrum around the FeKa complex energy (~6.4-7keV). 
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Fig. 4. The rms-spectra of ESO 113-GOlO, for light curve bin 
times of 5 ks (top) and 500 s (bottom). See text for details. 



This could be due to a reflected component remaining constant 
while the power-law continuum varies, as proposed for the two- 
component model (§3.3). An origin in a distant material is un- 
likely since no narrow 6.4 keV FeKc line is observed (cf. §|4|. 



0.3-1 keV flux (count s ) 

Fig. 5. 3-lOkeV versus 0.3-1 keV binned flux-flux plot of 
ESO 113-GOlO, compared with a best-fitting linear plus con- 
stant model (dashed line, representing two-component spectral 
variation) and a simple power-law (dotted line, representing high 
energy spectral pivoting). The two models are indistinguishable 
over the relatively narrow flux range covered by the data. 



3.3. Flux-flux analysis 

Many AGN show a characteristic steepening of their continuum 
X-ray spectra at higher fluxes. A useful insight into the origin 
of this continuum spectral variability can be obtained using flux- 
flux analysis. Tayl or et aP (|2003) showed how a simple plot of 
soft X-ray flux versus a harder X-ray flux can indicate whether 
spectral variability is due to pivoting at high energies (where 
the flux-flux plot takes a power-law form), or due to changes in 
the strength of a continuum component with a constant spectral 
shape plus a constant-flux, constant spectral shape harder com- 
ponent (where the flux-flux plot is linear with a positive intercept 
on the hard flux axis). Here we carry out a flux-flux analysis of 
the ESO 113-GOlO data, using the same combined MOS and PN 
events data used for the rms-spectral analysis. 

We first realize a flux-flux plot to compare the 0.3-1 keV and 
3-10 keV ranges. We measured the fluxes in simultaneo us 500 s 
bins, and following the approach of iTavlor et al.l ( l2003h . we bin 
the 3-10 keV flux as a function of the 0.3-1 keV flux, obtaining 
the standard error in the mean 3-10 keV flux of each bin from the 
spread in fluxes in the bin. The resulting flux-flux plot is shown 
in Fig.|5] A linear plus constant model provides a good fit to the 
data, with;^'^ =5.4 for 6 d.o.f., for positive constant off'set on the 
3-lOkeV flux axis of 0.082+0.014 counts"' (error bars are Icr). 
A single power-law (with no offset) also provides a good fit to 
the flux-flux plot (x^ =5.4 for 7 d.o.f.), for a power-law index of 
0.75 + 0.04. Thus, the existing data cannot distinguish between 
spectral pivoting of the continuum and a two-component origin 
for the spectral softening, primarily because the data span an in- 
sufficiently broad range of fluxes to see the characteristic change 
in gradient expected from the power-law. 

Although we cannot distinguish between spectral pivoting 
and a two-component spectral model, we can use the flux-flux 
method to determine the shape of any constant component un- 
der the assumption that th e two-component model is correct. 
Following the approach of iTavlor et al.l (l2003l) . we fit linear plus 



Energy (keV) 

Fig. 6. Fractional offset spectrum, showing fractional constant 
offsets (y-axis intercept divided by mean flux) for a linear plus 
constant model fitted to various energy bands versus the 0.3- 
0.5 keV reference band. In the case that the two-component 
model for spectral variability is correct, the offset spectrum 
shows lower limits on the strength of the constant component 
spectrum relative to the mean spectrum. 



constant models to the flux-flux plots measured from light curves 
obtained in a number of narrow bands. For the jc-axis, we use 
the 0.3-0.5 keV band as our 'reference' band. Thus all offsets 
measured in other bands are relative to this band only, under the 
assumption that there is zero flux in the constant component in 
the reference band. If there is a non-negligible flux in the ref- 
erence band, the implied offsets in other bands increase accord- 
ingly. Thus, if the two-component model is correct the offsets we 
measure from flux-flux plots here represent a lower-limit on the 
strength of any constant component. To remove the effects of 
the instrument spectral response, we normalize each measured 
offset by the mean flux in each band. The resulting spectrum of 
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Fig. 7. Spectra of ESO113-G010 (in the observer frame) observed in November 2005. A power-law has been fitted to the 2-5 keV 
data and extrapolated to lower and higher energies. A significant soft X-ray positive residual is seen below 1 keV, as well as the 
presence of a FeKo- complex near 6-7keV. For presentation, the data have been re-binned into groups of 15 bins for PN (10 for 
MOS), after group of a minimum of 100 and 50 counts per bin were used for the fit for the PN and MOS respectively. Left panel: 
0.2-12keV PN spectrum. Right panel: co-added MOS 0.3-lOkeV spectrum. 



fractional offsets is shown in Fig. |6] Since fractional offset in- 
creases significantly toward higher energies, the offset spectrum 
is harder than the mean spectrum, implying that any constant 
component is also significantly harder than the variable com- 
ponent, as expected if we are witnessing a constant reflection 
component. Notably the reflected flux contributes about 40% of 
the mean at the hardest energies, similar t o what is observed 
in MC G-6-30-15 using the same approach dVaughan & Fabiaiil 
12004 . 



4. Spectral analysis 

In this section, we present the spectral analysis of ESO 113- 
GOIO. In this article, we focus on the analysis of the 2-12 keV 
hard X-ray range (not affected by WA signatures), including the 
study of both the continuum and the Fe Ka line complex be- 
tween 6.4-7 keV. We co-added the MOSl and MOS2 data into 
a single spectral file to maximize the signal-to-noise ratio. The 
time-averaged PN and MOS spectra are binned to give a mini- 
mum of 100 and 50 counts per bin, respectively; while the sub- 
spectra used in the time-resolved spectral analysis (see below) 
are binned to give a minimum of 20 counts per bin. In all subse- 
quent fits (above 2 keV), we include the Galactic column den- 
sity (A^^^'=2.74 xlO^"cm"^), obtained from the coldens pro- 
gram using the compilations of Dickey & Lockman ( 1990). Note 
that all fit parameters are given in the source's rest frame, with 
values of //o=75kms"' Mpc"', and ^o=0.5 assumed through- 
out. The errors quoted correspond to 90% confidence ranges for 
one interesting parameter (A;i'^=2.71). Abundances are those of 
[Anders & Grevesse (1989). In the following, we use the updated 
cross-sections for X-ra y absorption by the inters t ellar m edium 
(tbabs in XSPEC) from lWilms. AUen. & McCravl (|2000|) . 

4. 1 . The main spectral components 

As a first step, we fit an absorbed power-law model over the 
2-5 keV energy range where the spectrum should be relatively 
unaffected by the presence of a broad soft excess, a warm 
absorber-emitter medium, an Fe Ka emission line, and a con- 
tribution above 8keV from a high energy Compton reflection 



hump. In this energy range, the PN data are well fitted by a single 
power-law model with r= 2.02+0.04 Ct-^/d.o.f. =228.3/214). For 
comparison the value from the 2001 observation is 1.96±0.22. 
This power-law index is consistent with values found for broad- 
line radio-quiet quasars: r=1.90 with a standard deviation of 
0.27 (Porquet et al. 2004a). The 2-lOkeV time-averaged lumi- 
nosity (inferred from the parameter fit reported in Table [T]i is 
4. Ox 10"^^ erg s"' . Figure|7]displays the 2-5 keV PN spectrum ex- 
trapolated over the 0.2-12keV broad band energy. A strong pos- 
itive residual is seen below 1 keV due to the presence of a soft 
X-ray excess; and in addition, there is a positive residual in the 
6-7 keV energy range (see § I4.2l i. For comparison, we also re- 
port the co-added MOS 1 and MOS 2 spectrum, which confirms 
the presence of a complex soft excess as well as a positive devia- 
tion near the 6-7 keV. We henceforth focus only on the PN data, 
since the PN CCD has a better sensitivity over a broader energy 
range (0.2-12 keV) compared to the MOS CCD. However we 
have checked that all the spectral fits are consistent with those 
obtained with the MOS, albeit with much lower photon statis- 
tics. We would like to note that the RGS data show the presence 
of several narrow blue-shifted absorption lines that reveals the 
presence of a warm absorber (WA) in outflow. The presence of 
this medium has to be taken into account in order to fit properly 
the overall 0.2-12 keV PN spectrum. Spectral analysis combin- 
ing RGS and PN data will be presented in a forthcoming paper. 

4.2. The Fe Ka line complex 
4.2.1. Time-averaged analysis 

The positive residuals between 6.5-7 keV in the quasar frame 
can be most probably a ssociated with a FeKg complex as 
observed in most AGN (Porqu et et al.l l2004at iPiconcelli et al.l 
1.2005; Jimenez-Bailon et al. 

2005r The mean iron line profile 
as a ratio of the continuum above 2keV is displayed in Fig. [8] 
and shows a double-peaked line profile. Therefore, we fit the 2- 
12 keV energy range with an absorbed power-law continuum and 
two narrow Gaussian line profiles (Table [TJ. We found a much 
better fit, with Ax^=33.2 for 4 additional parameters (F-test > 
99.9%). The fine energies are 6.50!°°^ keV and 7.00+0.04 keV. 
The F-test probabilities for each line are 99.7% and 99.98% for 
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Table 1. Best-fitting spectral parameters of PN time-averaged spectrum in the 2-12 keV energy range with an absorbed (Galactic, 
A^H=2.74xlO^**cm"^) power-law (PL) component plus a line profile: zgauss: Gaussian profile; and diskline: profile line emitted 
by a relativistic accretion disc for a non-rotating black hole (Fabian et al. 1989) assuming an emissivity law q equal to -2. (a): 
^out=l 000 7?g. (b): /?out=l-2 X Ria. The line fluxes are expressed in units of 10 ''phcm^- s"'. F-test probabilities are calculated with 
the corresponding power-law model as reference. Fixed parameters are indicated by (f). 



Model 


r 




Line parameters 




X^ld.o.f. 


F-test 






E 




e 


Flux 


EW 








(keV) 


(eV) (^g) 


(deg) 




(eV) 





PL 

PL + 2 X zgauss 
PL + diskline'"' 

PL + ANNULUS**' 



1.95+0.02 
1.97±0.02 

2.00+0.02 
1.97±0.02 



6 50+°"'' 
7.00±0.04 
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Fig. 8. Fe K line profile (quasar frame) when PN time-averaged 
data are fit with a power-law continuum. For presentation, the 
data have been re-binned into groups of 10 bins for PN, after 
group of a minimum of 100 counts per bin were used for the fit. 
The grey vertical line corresponds to 6.4 keV. 



the 6.50 keV and 7.0 keV lines, respectively. Allowing the line 
widths to vary does not significantly increase the goodness of 
fit, with - 0.5 for two additional parameters. 

The former energy is compatible within the error bars with 
a line emitted by moderately ionized iron (~Fexix-Fexxii). 
However this line is more likely a blue or red peak of a rel- 
ativistic line, indeed a genuine Fexix-Fexxii (L-ions) emis- 
sion line complex would sufl'er fro m resonant Auger destruc- 
tion dRoss et al.lll996tlLiedahill2005l) . The latter line energy cor- 
responds to highly ionized iron line (Fexxvi) or to a highly 
blue-shifted iron line from lower ionization matter However 
there is no evidence for narrow emission from neutral iron, the 
90% confidence upper- limit to a narrow (cr=10eV) iron line at 
6.4 keV (a 6.5 keV narrow line with a width of 10 eV is also fit- 
ted to take into account its emission contribution) is only 32 eV. 
If the 6.5 keV line contribution is not taken into account, we 
would infer £'W(6.4keV) <47eV. From the branching ratio of 
Fe_K6/Fe Ko-, expected to be less than 0.145 (Palme rTet al.l 
I2003L and references therein), the maximum contribution of the 
Fe K/S (from cold iron) to the 7.00 keV line is only up to 5 eV. 



We also tested if these two features can be explained by a 
double-peaked relativistic line profile emitted by a r elativistic ac- 
creti on disc for a non-rotating black hole (diskline; iFabian et alj 
Il989i) . The improvement of the fit is also significant (F-test > 
99.9%) and the line energy is about 6.8 keV, i.e., corresponding 
to a highly ionized iron ion (Fe xxv-Fe xxvi). Then, we checked 
whether the double-peaked like profile can be explained by emis- 
sion from a disc annulus (using the diskline profile with the 
width fixed to 20% of Ri^, i.e. Roui-l-2 x Ri^)- We fixed the 
disc inclination to 45° (representative of an intermediate type 
Seyfert). We found a very good fit to the data with a line energy 
of 6.79+0.05 keV and an inner annulus radius of 191+^3 7?g, in- 
dicating that the line is not emitted in the inner part of the disc. 
In the case where the disc inclination is not fixed to 45 degrees, 
then the line solution is degenerate, whereby the inner annulus 
radius increases with the disc inclination, while the line energy 
decreases. 

Monte Carlo simulations were done to assess the significance 
of detecting the emission complex in this time-averaged spec- 
trum. These simulations, described in detail in the Appendix, 
demonstrate that regardless of the model used, the emission 
complex is indeed significantly detected. 

We also checked for an upper limit to a line at 5.39 keV, since 
such a line was found in the 2001 XMM-Newton observation 
of May 2001. We fixed the energy of the line to 5.39 keV (see 
below) , and the line width to 0. 1 keV (as done in iPorquet et alJ 
'2004t^. We find EW<22eV (at 90% confidence level for one 
interesting parameter). 

For comparison, we have re-processed the XMM-Newton 
May 2001 observation with the same version of the SAS 
(SAS 6.5.0) used for the 2005 observation. We fit the data 
with an absorbed (Galactic column density) power-law in 
the 2-lOkeV energy range and four Gaussian lines in order 
to take into account the emission contribution of the lines 
near/at 5.4 keV, 6.4keV, 6.5 keV, and 7.0keV. The three lat- 
ter line energies are fixed as well as their widths to 10 eV. 
Only the energy of the first Gaussian line is le t free t o vary 
with t he width fixed to 0. 1 keV (as done in iPorguet et al.l 
l2004bl) . One should notice that in iPorquet et aH (l2004bl) . 
only the line at 5.4 keV was fitted. We found for the four 
lines: £Wi=246±147 eV (£,=5. 39+0.09 keV), EW2 <190eV 
(£2=6.4keV), EW3 <100eV (£2=6.5 keV), and EWa <208eV 
(£2=7.0keV). The EW for the line at 6.4keV, 6.5 keV, and 
7.0 keV are compatible with the values found during the 2005 
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Fig. 9. Energy-time map. The energy values are in the quasar 
frame. The contours represent 1.5^.0 sigma above the contin- 
uum (in 0.5 increments), with darker bins representing an excess 
above the continuum. The segments 1-3 (delimited by red ver- 
tical lines) correspond to the time selections used for the time- 
resolved spectroscopic analysis. See text for explanation. 



observation. For the line at ~5.39keV, the EW are not compat- 
ible and would mean that if the line detected during the 2001 
observation is real, it is transient. 



4.2.2. Time-resolved analysis 

To test for possible rapid shifts in the energy and/or flux of 
the Fe Ka emission complex, we have created X-ray intensity 
maps in the energy-time plane using the PN data. Photons from 
the source cell were accumulated in pixels in the energy-time 
plane. The pixel distribution was smoothed in energy by the 
instrumental resolution, using a Gaussian of 140 eV (appropri- 
ate for the single plus double events with the latest calibra- 
tion), and smoothed in time using a top-hat function of width 
20 ks. Each time-slice was background-corrected by subtract- 
ing a time-dependent background spectrum measured in an off- 
source region on the same detector chip as the source. The source 
continuum was modeled as an absorbed power-law, of variable 
amplitude and slope but time-invariant absorption column den- 
sity. This continuum was subtracted, leaving positive and nega- 
tive residuals that comprise noise plus any emission or absorp- 
tion components on top of t he continuum. Mor e information 
about this method is given in iTurner et al.l d2006h . The "signal- 
to-noise" (S/N) map presented in Fig.|9]is the ratio of the fluctu- 
ation amplitude to the calculated noise. 

The highest S/N observed is 4 and there are several points of 
significant excess emission above the continuum. However, we 



5 6 7 8 

Rest Energy (keV) 

Fig. 10. The PN data/model sub-spectra of ESO113-G010 
(quasar frame) for the three time intervals. The data have 
been fitted assuming a power-law underlying continuum. Upper 
panel: 0-40 ks sub-spectrum (segment 1). Middle panel: 40- 
70 ks sub-spectrum (segment 2). Lower panel: 70-100 ks sub- 
spectrum (segment 3). 



note here that the energy-time map highlights fluctuations above 
the continuum which may be real spectral features, but does not 
assess the significance of variations in those features, assuming 
they are real. The apparent fluctuations in the energy-time map 
could be consistent with the complex emission features detected 
in the time-averaged spectrum being constant, e.g. if the average 
significance in one independent 10 ksec segment is 2 sigma 
above the continuum, then we might expect statistical variations 
from 0^ sigma in a 100 ksec observation. 

To quantify whether the line features are variable, we split 
the data set into three time intervals, 0^0 ksec, 40-70 ksec 
and 70-100 ksec, which are chosen to maximize the diflFerences 
in the Fe K line region which are highlighted by the energy- 
time map. Note that in the case that the lines are constant, this 
approach essentially maximizes the eff'ect of noise fluctuations 
that are picked up by the energy-time map, so that the statis- 
tical tests we carry out are in some sense the most optimistic 
for detecting line variability. The line profiles for each segment 
(data/model ratio for a power-law continuum over the 2-lOkeV 
energy range) are shown in Fig. [10] We then fit these three sub- 
spectra with (individual or double) Gaussian line profiles; the re- 
sults are reported in Table |2] Though the line profiles in Fig.fTOl 
appear to be diff'erent, adding a double Gaussian line profiles at 
6.5 keV and 7.0 keV (as found in the time-average spectrum) are 
significant (F-test significance=99. 5-99.97%) in all segments. 
We also carried out Monte Carlo simulations to verify the sig- 
nificance of lines fitted to the individual sub-spectra; we detail 
these results in the Appendix (see IA.3l l, but note here that they 
confirm that the lines are significantly detected in the individual 
sub-spectra. 

Having verified the statistical significance of the lines in each 
of the 3 segments, we check for significant variability of the line 
flux. First, we fit simultaneously all segments with a power-law 
model plus two Gaussian lines with line energies fixed to 6.5 keV 
and 7.0 keV and line widths and line normalisations are free to 
vary within each segment. The contour plots at 90% confidence 
level for two interesting parameters (A;^'^=4.61) are reported for 
each line in Fig.[TT] and show that assuming two Gaussian lines 
there are no statistically significant line variations. We also sim- 
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Table 2. Best-fitting spectral parameters of PN time-resolved spectra in the 2-10 keV energy range with an absorbed 
(Galactic, A^H=2.74xlO^''cm"^) power-law (PL) component plus a Gaussian profile (zgauss): . The line fluxes are expressed in 
10"^ ph cm"^ s~^ . F-test probabilities are calculated with the corresponding power-law model as reference. (*) The F-test probabil- 
ity for freeing the line width in the fit is 92%. (a) Monte-Carlo simulations have been performed (see text). Fixed parameters are 
indicated by (f). 
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Fig. 11. Width and normaUsation line contour plots at 90% confidence level (A;^f^=4.61) for segment 1 (dotted-Une), segment 2 
(continuous line), and segment 3 (dashed-line). Left panel: line at 6.5 keV. For segment 2, only an upper limit is found, corresponding 
to the continuous line. Right panel: line at 7.0 keV. For segments 1 and 3, only an upper limit is found. 



ulated 100 spectra using the fakeit command, and assuming an 
initial model consisting of a power-law plus a single Gaussian 
emission line at either 6.5 or 7.0 keV. We assumed that the line is 
intrinsicaUy constant, with parameters described by values, and 
used an exposure time of 29.8 ks. We found that the observed 



values of line intensity and line width cr were usually consistent 
with the resulting spread in 100 simulated values. This suggested 
that we could not rule out at high significance the hypothesis that 
these parameters are intrinsicaUy constant. 
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5. Discussion 

We report the data analysis results of a long lOOks XMM- 
Newton observation in November 2005. We have used 
both model-independent techniques (e.g., rms-spectra, flux-flux 
plots) and spectral analysis (time-averaged and time-resolved 
spectroscopy). 

5.1. Timing analysis 

Light curves reveal that this source was highly variable on short 
timescales (~500s), and was on average in a higher/softer state 
compared to the May 2001 observation. The Power Spectral 
Density (PSD) analysis indicates (at 95% confidence) a break 
at a frequency of 3.7(^[ ") x 10""* Hz. This measured cut-off" fre- 
quency is comparable to those measured in some other rapidly- 
variable Seyferts, such as MCG-6-30-15 and NGC4051. For a 
range in L[,oi from lO'*'*'^ - 10'*^ erg s"' (i.e. assuming a min- 
imum equal to th e sum of FIR, op tical, and X-rays luminosi- 
ties, reported by Pietsch et al.ll 19981 and in this work), we infer 
from the mass-luminosity-time-scale relation of iMcHardv et al.l 
( 1200 6) that Mbh ranges from 4 X 10*'-10^ Mq and the source is 
accreting at or close to the Eddington rate (assuming even higher 
Lboi implies that the AGN is a super-Eddington accretor). 

The rms-spectra obtained using light curves binned to 5 ks 
and 500 s show that Fvar appears to decrease above about 2 keV 
for both bin timescales. Therefore the soft X-ray excess down 
to 0.3 keV appears more variable compared to the hard en- 
ergy range. This is different from the typical rms-spectra from 
other Seyferts, which tend to show a broad peak in the vari- 
ability amplitude between 0.8 and 2keV, and a much lower 
variabihty amplitude below 0.8 keV (e.g., iFabian etld] 120021: 
llnoue & Matsumotoll2003l:lPonti et aLll2006 l). From the flux-flux 
analysis, the existing data cannot distinguish between spectral 
pivoting of the continuum and a two-component origin for the 
spectral softening, primarily because the data do not span a 
broad enough flux range to see the characteristic change in gra- 
dient expected from the power-law. However, in the case of the 
two-component model being correct for this source, the frac- 
tional off'sets measured in the flux-flux plots increase signifi- 
cantly toward higher energies (similarly to what is observed in 
MCG-6-30-15), as expected if we are witnessing a constant re- 
flection component. Because there is no evidence for distant, 
cold material such as a torus (§4), then if the Compton reflec- 
tion component is real, the reflection is more likely associated 
with the accretion disc. 

5.2. Spectral analysis: The FeKa line complex 

Instead of a single highly redshifted line at about 5.4 keV 
as observed in the 4ks May 2001 observation, we found in 
the November 2005 time-averaged observation two narrow 
FeKa emission lines at about 6.5 keV and 7.0 keV. The former 
rest-frame energy is compatible within the error bars with a line 
emitted by moderately ionized iron (~Fexix-Fexxn). However 
this line is more likely a blue or red peak of a relativistic line, 
indeed a genuine Fe xix-Fe xxii (L-ions) emission line complex 
would suffer from resonant Auger destruction dRoss et aLl 
[19961 iLiedahll |2005|) . The latter line energy corresponds to 
highly ionized iron line (Fexxvi) or to a highly blue-shifted 
iron line from lower ionization matter. The two features can 
be well represented by a double-peaked line profile from the 
disk with a rest-frame energy of 6.75-6.80keV. In case of 
an annulus diskline profile an inner radius of about 200 7?g is 



inferred. Alternatively the Fexxvi line could be the signature 
of a strongly photoionized (log^ >3), circumnuclear matter 
(seen in emission) with an iron overabundance of a factor 3, 
and/or very high colum n densities with A^h warm >> 10^^ cm"^ 
(iBianchi & Mattll2002l) . There is no significant evidence for line 
variability, however with the current data we also cannot reject 
the possibility that the lines do vary in response to continuum 
variations. Further XMM-Newton monitoring of this source can 
help to establish or reject line variability at different time-scales 
(from day(s) to months) and also allow us to determine if 
the line profile is indeed double-peaked and originates in the 
disk or if there exist independent or transient red-/blue-shifted 
emission lines. Shorter time-scales will only be reached by 
higher collecting area (and resolution at the iron Kor energy 
band) instruments such as Constellation-X and XEUS. 

In the May 2001 observation of ESO 113-GOlO we observed 
a line feature at 5.39 keV, which was detected at 99% confidence 
(based on Monte Carlo simulations) with EW= 233 + 140 eV. 
We notice that a weak feature is seen at 5.4 keV near 20 ks in the 
energy-time map (Fig. 9), but this line has a detection probabil- 
ity of only 38% using Monte-Carlo simulations. The 90% upper 
limit on EW for a 5.39 keV line in the total observation is 22 eV 
(Section 4.2.1). Therefore, if the line observed in May 2001 is 
real, it is transient. We note however that, if we treat the total data 
set on ESO 1 13-GOlO as a set of independent ~ 4 ksec segments, 
the chance of seeing a 99% ffuctuation in a single segment by 
chance becomes quite high. However, in support of the possibil- 
ity the 5.39 keV feature was real, we also point out that the ffux 
and spectral hardness in May 2001 corresponds to the extreme 
low and high values observed in 2005 (see Fig. 1), which might 
be more favourable to transient line formation (though we also 
note that such a line feature was not observed during the 2005 
ffux minimum). 

In addition, as suggested by the 2001 observation, no sig- 
nificant constant narrow 6.4 keV FeKa line is observed with 
an upper-limit of < 32 eV, hence excluding any dominant emis- 
sion from distant cold matter such as a torus. For a neutral so- 
lar abundance reffector subtending 2:7rsr"' to the line of sight, 
illuminated by a F = 1.9 continuum and inclined at 45°, the ex- 
pected iron line equi valent width is ~ 125 eV with respect to the 
primary continuum (lGeorge&Fabianlll99lh . Thus in ESO 113- 
GOIO any cold, Compton-thickreprocessor outside of the line of 
sight is likely to subtend a much smaller solid angle, of the order 
< 7r/2sr"'. 

The optical classification as a Sey fert type 1 .8 of ESO 113- 
GOIO (observed in November 1996: iPietsch et all [T998h sug- 
gests the presence of a reddened broad line region (BLR). 
For a reddening of E(B-y) = 0.8-1.0, typical for a Seyfert 
1.8 (IVeron-Cetty & Veronll2006l) . the hydrogen column density 
would be equal to or greater than about 4-6x10^' cm"^ (as- 
suming a dust-to-gas ratio equal to or smaller than the Galactic 
value; iBohlin et al.l Il978l) . However, during the present 2005 
XMM-Newton observation (also suggested by the 2001 obser- 
vation), instead of an intrinsic neutral absorption (in addition to 
the Galactic column density value of 2.74 x 10^** cm"^) expected 
about below 1 keV, a soft excess is observed. A resolution of 
this apparent discrepancy, without the need to invoke a very high 
dust-to-gas ratio is either that dust that obscures the BLR clouds 
could lie out of the direct line of sight toward the X-ray emit- 
ter; or dust responsible for the BLR reddening could be mixed 
in with the warm absorber-emitter medium detected in the RGS 
spectrum. Alternatively, the apparent disagreement between the 
(non-simultaneous) optical and X-ray classifications could be 
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explained by a time variability of th e optical and/or X-r ay ab- 
sorbing properties (e.g., H1320-H551: lBarcons et al.ll2003h . 
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Appendix A: Monte-Carlo simulations 

A.1. Method 

We carried out rigorous tests of the significance of the lines in the 
time-averaged spectrum and in each segment using Mont e Carlo 
simulations; see Porquet etal. (2004b) and Markowit z et al.l 
(l2006h for full details. For the null hypothesis, we assumed that 
the spectrum is simply an absorbed power-law continuum, with 
the same parameters as the absorbed power-law model fitted to 
the real data in each segment. We used the XSPEC FAKEIT 
command to create 1000 fake EPlC-pn spectra corresponding to 
this model, with photon statistics appropriate for each exposure 
(38.6, 24.0, and 26.9 ks for segments 1, 2 and 3, respectively). 



and grouped each spectrum to a maximum of 20 counts per bin. 
For each faked spectrum, we re-fit the null hypothesis model, 
yielding a "modified" null hypothesis model. We ran FAKEIT a 
second time using this re-fit model, again with the appropriate 
exposure time; this process accounts for the uncertainty in the 
null hypothesis model itself. For each new spectrum, we re-fit 
the null hypothesis model, recorded an initial value, and then 
added a Gaussian component to the fit, as described below. The 
line centroid energy was constrained to the 4.5-8.0 keV range, 
and the line normalization was allowed to be positive or nega- 
tive. We stepped the Gaussian centroid energy over this range in 
increments of 0. 1 keV, fitting separately each time to ensure the 
lowest value was found. For each spectrum, we compared the 
lowest value with that from the null hypothesis fit, to obtain 
1000 simulated values of the A^', which we used to construct 
a cumulative frequency distribution of the Ax^ expected for a 
blind line search in the 4.5-8.0 keV range, assuming the null 
hypothesis of a simple power-law with no line is correct. 

A.2. Results for the time-averaged spectrum 

The observed 6.5 and 7.0 keV lines had Ax^ values of 13.2 and 
19.0, significant at 97.4% and 99.7% confidence, respectively, 
according to the simulations. However, since we have detected 
two lines that are close in energy, we can quantify their detec- 
tion significance as a pair The likelihood that both lines are due 
to photon noise is (l-0.974)( 1-0.997) = 8x10"^ i.e., the pair 
is significant at just over 99.99% confidence. As an alternative 
quantification of the detection significance of the pair, we ran 
simulations which tested for the presence of a pair of narrow 
Gaussians with a fixed energy separation of 0.50 keV and with 
both line intensities left free. The observed Ax^ value of 33.2 
is significant at >99.9% confidence. However, as there was no 
a priori expectation for the value of the energy separation, us- 
ing such a line profile tests only a small subset of the full range 
of line pair profiles possible, and so this probability should be 
treated as an upper limit. Assuming the annulus profile is cor- 
rect, the observed Ax^ value of 28.8 corresponds to a detection 
probability of >99.9%. 

A.3. Results for the time-resolved sub-spectra 

Segment 1 (0-40ks): The observed 6.51 keV line had a Ax^ of 
9.4, significant at 85.8% according to the Monte Carlo simu- 
lations. The observed 7.03 keV line has a Ax^ of 4.2, signifi- 
cant at only 18.6% confidence, i.e., we cannot significantly rule 
out the hypothesis that the 7.03 keV feature is due to photon 
noise. However, the likelihood that both lines are due to photon 
noise, calculated in the same manner as for the time-averaged 
spectrum, is 0.116, i.e., the pair of lines is significant at 89.4%. 
Performing Monte Carlo simulations assuming two lines with a 
fixed energy separation yields (an upper limit of) 92.0% confi- 
dence. 

Segment 2 (40-70 ks): The simulations suggest the broad 
Gaussian 7.0 keV line, with an observed Ax^ of 17.9, is signifi- 
cant at 99.7% confidence. 

Segment 3 (70-lOOks): The simulations suggest the broad 
Gaussian 6.7keV line, with an observed Ax^ of 13.9, for only 
3 additional parameters is significant at 99.9% confidence. 
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